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We study the local correlations in the super Tonks-Girardeau gas, a highly excited, strongly 
correlated state obtained in quasi one-dimensional Bose gases by tuning the scattering length to large 
negative values using a confinement-induced resonance. Exploiting a connection with a relativistic 
field theory, we obtain results for the two-body and three-body local correlators at zero and finite 
temperature. At zero temperature our result for the three-body correlator agrees with the extension 
of the results of Cheianov et al. [Phys. Rev. A 73, 051604(R) (2006)], obtained for the ground-state of 
the repulsive Lieb-Liniger gas, to the super Tonks-Girardeau state. At finite temperature we obtain 
that the three-body correlator has a weak dependence on the temperature up to the degeneracy 
temperature To- We also find that for temperatures larger than Tu the values of the three-body 
correlator for the super Tonks-Girardeau gas and the corresponding repulsive Lieb-Liniger gas are 
rather similar even for relatively small couplings. 

PACS numbers: 



I. INTRODUCTION 

The experimental capability of tailoring tightly confined trapping potentials opened the way to realizing 
one-dimensional (ID) interacting systems with ultracold atoms [1,2]: when the transverse motion of atoms 
is confined to zero-point oscillations, the effective Hamiltonian describing their equilibrium properties and 
their dynamics is ID. It is then possible to experimentally simulate paradigmatic ID many-body models: 
ID interacting Fermi gases provide an experimental realization of the Gaudin-Yang model [3], while ID 
Bose gases are very well described by the Lieb-Liniger (LL) model (see the reviews [4, 5]). These models 
have been theoretically studied for decades since they display a rich variety of non mean-field features: the 
study of their equilibrium properties motivated the developments of analytical and numerical techniques 
and, at the same time, they were used as a benchmark for testing non-perturbative techniques. 

The experimental realization of the LL model with ultracold bosons not only renewed the interest 
in the equilibrium and dynamical properties of the model, but also called for the study of its excited, 
strongly correlated states. Indeed, tuning the effective ID scattering length it is possible to vary the 
coupling constant 7 of the LL model: 7 can be made very large close to a confinement-induced resonance 
[6]. For large positive 7 one approaches the Tonks-Girardeau (TG) limit [7], while suddenly switching 
to the other side of the resonance one can prepare a highly excited many-body state, the super Tonks- 
Girardeau (STG) state [8]. Exploiting a confinement-induced resonance and using a gas of Cesium atoms, 
the STG state has been recently realized [9]: the crossover from the TG to the STG regime has been 
studied by determining the collective mode frequencies and the dynamics through the crossover has been 
characterized by measuring the particle loss and the expansion [9]. 

Several properties of the STG gas have been discussed recently in the literature [8-12]: in [10] it was 
shown that the STG gas corresponds to a highly excited state in the Bethe ansatz solution of the LL model 
with attractive interactions (i.e., negative 7), characterized by real Bethe roots. In [11] the connection 
between the STG gas and a, ID hard sphere Bose gas (with hard sphere diameter almost equal to the ID 
scattering length) was discussed. The realization of effective STG gases in strongly attractive ID Fermi 
gases has also been proposed and investigated [13, 14]. The dynamics in the crossover between the TG 
and the STG gases when the coupling constant 7 is suddenly quenched from a positive value (in which 
the system is in the ground state) to a negative value was considered in [12, 15]: after the quench the 
system is in a metastable state, however, solving the exact dynamics it is possible to see that such a state 
can be stable for rather long times (~ 100ms). The metastability of the STG gas can be understood in 
terms of the very small overlap between the ground state of the TG gas and the collapsed cluster states. 

In this paper we focus on the computation of local two- and three-body correlations of the STG gas, 
as these important quantities determine the rates of inelastic processes, such as photoassociation in pair 
collisions and three-body recombination. We use a recently introduced method [16] which allows for 
the determination of equilibrium expectation values in the repulsive LL gas: extending this approach 
to the STG state, we can compute two- and three-body correlations in such a highly excited, strongly 
correlated state both at zero and finite temperature. At T = we find that our results for the three-body 
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correlator g 3 are in agreement with the extension to the STG state of the results of [17], obtained for 
the ground-state of the repulsive Licb-Liniger gas. At finite temperature, g% displays for intermediate 
values of the coupling constant a weak dependence on the temperature up to the degeneracy temperature 
Tb; furthermore for T > Tb the ratio between the three-body correlation of the STG gas for a coupling 
constant 7 (with 7 < 0) and the corresponding value computed at equilibrium for the LL gas with 
coupling constant |7| is rather close to one even for relatively small values of I7I (i.e. for I7I > 20). 

The plan of the paper is the following: in Section II we review basic facts of the LL model and introduce 
the STG state, while Section III is devoted to the computation of local expectation values in the STG 
gas. In Section IV we show our results for zero temperature, and the findings for finite temperature are 
discussed in Section V. In Section VI we present our conclusions. 



II. THE SUPER TONKS— GIRARDEAU GAS 



The LL Hamiltonian describes N non-relativistic bosons of mass m in one dimension, interacting via 
a two-body (^-potential [18]: 



h 2 N d 2 

* = -2^Eft£+2Aj>*-*i)- (1) 

The coupling A of the II? 5-function contact potential in the Hamiltonian (1) can be expressed in terms 
of the parameters of the three-dimensional Bose gas [6] as 

A = ^ * . , (2) 



where a^u is the three-dimensional scattering length of the Bose gas [1, 2], a± = y/ h/mojj_ is the harmonic 
oscillator length of the transverse confinement with trap frequency wj_ and C ~ 1.0326 is a constant. The 
effective coupling constant of the LL model is given by the dimcnsionless quantity 

2mA 

7 = -PT~ . 3 
nrn 

where n = N/L is the density of the gas (L is the length of the system). 

The limit 7 -C 1 is the weak coupling limit where the Bogoliubov approximation gives a good estimate 
of the ground-state energy of the system [18] . For large positive 7 one approaches the TG limit [7] , where 
the combined effect of the reduced dimensionality and the strong repulsion leads to an effective Pauli 
exclusion and the ground-state wave function can be mapped to the wave function of free fcrmions. A 
value of 7 w 5 was reached in [19], while in [20], using an additional shallow optical lattice along the 
longitudinal direction, effective values sa 200 for 7 were achieved. As one can see from (2), by tuning 
a-3D ~ a± one can have large 7 and pass from positive to negative values of 7 [6, 9]. 

In the LL model temperatures are usually expressed in units of the quantum degeneracy temperature 

h 2 n 2 , s 

in the following we use the scaled temperature 

T 



Notice that defining the thermal De Broglic wavelength At = \J2ttH 2 jvnk^T [1, 2], one has tiAt d = 
showing that for T ~ degeneracy effects arise. For ID tubes having N ~ 100 atoms and size L ~ 10/^m 
(corresponding to longitudinal frequencies uj z ~ 2ix ■ (1 — 5) Hz), the degeneracy temperature Td is 
~ 300nK. For a transverse trapping confinement uj± ~ 2ir ■ 5kHz, one has Hu>±/2 ~ fee ■ lOOnK, and then 
scaled temperatures as low asT~ 0.3 are realistically reachable. 

As shown by Lieb and Liniger in their original paper [18], the eigenvalue problem of the Hamiltonian 
(1) can be solved in terms of a coordinate Bcthc ansatz. The equations of motion are just free Schrodinger 
equations in the domain where the coordinates of the particles are all distinct. If we denote by R\ the 
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subset of the configuration space where X\ < x 2 < ■ ■ ■ < xjy , the solution of the equations in Ri is given 
by the Bethe wave function 

XN (xi,x 2 ,...,x N ) = Y,a(P)e i ^ p ^ > ( 6 ) 
p 

where ^ p denotes a sum over permutations of the wavevectors . . . , k n } characterizing the state. 
For configurations outside Ri the solution is easily obtained using the symmetry of xn with respect to 
the Xi. For the permutations P : (fe, I, fc Q3 , . . . , k aN ) and Q : (Z, k, fc Q3 , . . . , k aN ) the relation between the 
coefficients in the sum appearing in (6) is given by 

a(Q)= k ~ l ~ iK a(P), (7) 
k — I + in 

where k = Tpr A = wy. Hence the wave function gets multiplied by the factor a(Q) / a(P) whenever two 
particles with momenta p\ = hk and p 2 = hi are exchanged, therefore the two-body S'-matrix of the 
Lieb-Liniger model is 

S^{k,X)= l —^ 1 (8) 
k + ijn 

where hk = hk\ — hk 2 is the momentum difference. 

With periodic boundary conditions the momenta in (6) are constrained by the Bethe equations 

^ L --n !rt i+ — i=i,.~,N, (9) 

f-J- kj - ki - i-yn 

and the ground-state energy is given by E = ^ S^Li kj- The roots of these equations for the ground 
state are all real for 7 > [21] and the energy for large coupling 7 ^> 1 is given by 

E TT 2 h 2 n 2 ( 2Y 2 
N 6m \ 7 / 

The thcrmodynamical properties of the LL model for 7 > can be obtained by the thermodynamic 
Bethe ansatz (TBA) , as shown originally by Yang and Yang [22] : in Appendix A we summarize the TB A 
equations for the repulsive LL gas. 

For negative coupling, 7 < 0, the ground-state is a cluster-like state having complex Bethe roots with 
energy [23, 24] 

E = -^-N{N 2 ~l). (11) 

Due to the attraction, the particles tend to collapse in the same region of the space and the system 
does not have a well-defined thermodynamic limit: the ground-state energy per particle E/N — > —00 for 
N -)• 00. 

However, the existence of a stable gas- like state has been proposed [8]: the STG state. This state is 
an eigenstate of the attractive (7 < 0) LL model characterized by Bethe roots which are all real. It 
was shown [10] that this is a highly excited state in the attractive regime which is, however, stable in a 
wide range of coupling strength. As it was recently confirmed in the experiment [9], the STG state can 
be created from a TG gas by an abrupt change of the sign of the interaction. The large kinetic energy 
inherited from the TG acts like a Fermi pressure and it cannot be quenched instantly. Moreover, for large 
coupling strengths the two wave functions are almost identical and their overlap is very large. A detailed 
analysis of the root distribution can be found in [10, 12]: it is found that the configuration of the real 
roots corresponds to the roots of the repulsive LL gas but with the sign of the coupling changed. 

These results suggest a simple way to calculate the energy and other quantities for the STG gas: one 
can use the equations and formulae for the repulsive (7 > 0) LL gas for 7 < 0. In this way in [8, 10] the 
following asymptotic expression for the energy of the STG gas with coupling constant — 17| (with I7I ^> 1) 
was obtained: 

e(7< 0)«£(l- A)" 2 , (12) 
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FIG. 1: Energy eipf) for the STG (upper lines) and LL (lower lines) gas. The solid lines are the exact values, while 
the dotted lines refer to the asymptotic expressions (10) and (12), the dotted-dashed line shows the asymptotic 
value 7r 2 /3. 

where e(j < 0) = E(^y < 0)/NksTD [sec cqn (A7)]. The corresponding asymptotic expression 
for the ground-state energy of a repulsive LL gas with coupling constant 7 > is (for 7 3> 1) 

e (l) ~ if + 7j ) j as can be seen directly from eqn (10). 

The energy of the STG gas at finite values of — I7I can be obtained at T = using the integral equations 
reported in Appendix A: in these equations, valid for the ground-state of the LL model with 7 > 0, one 
has to replace 7 with — ItI * The result is plotted in Fig. 1, where we show for comparison the exact 
energies calculated from the LL integral equations for the STG and LL gases. 

To simplify the notations in this figure as well as in the following figures and sections, for the STG 
results the coupling 7 should be understood as the positive parameter —7 = I7I, while the corresponding 
LL results refer to equilibrium values for the LL gas with 7 > 0. 

III. LOCAL CORRELATORS IN THE SUPER TONKS— GIRARDEAU GAS 

In this section we present a discussion of the way to compute local correlators in the STG gas both at 
zero and finite temperatures. The quantities we are interested in are defined as 



where k = 1, 2, • • • . Equilibrium local correlators for the repulsive LL gas were discussed in several papers: 
exact results based on the Yang- Yang equations and the Hellmann-Feynman theorem are available for 
two-body correlations [25, 26], while the local three-body correlation g^ was determined at zero temper- 
ature in [17]. Asymptotic expressions for for large and small 7 were presented in [25]. Results for c/3 
at finite temperature were obtained in [16]. 

To compute the correlators gu in the STG we use the method recently introduced in [16]: this approach 
is based upon the fact that in (1+1) dimensions the repulsive LL model can be obtained as a suitable 
non-relativistic limit of an integrablc relativistic quantum field theory, the sinh-Gordon (sh-G) model. 
This model is defined by the Lagrangian density 



where = (f>(x,t) is a real scalar field, c is the speed of light and the parameter ll is related to the 
physical mass m by li 2 = iram 2 c 2 /h 2 sin(7ra), where a = hcg 2 /(8ir + heg 2 ) [27]. The energy E and the 
momentum P of a particle can be written as E = mc 2 cosh P = mcsinh^, where 6 is the rapidity. 
Since the sh-G dynamics is ruled by an infinite number of conservation laws, all its scattering processes 
are purely elastic and can be factorized in terms of the two-body S-matrix [27] 
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FIG. 2: Difference between our evaluation of 31(7) and the exact value g\ = 1 as a function of 7 at T = for 
the STG gas. The solid lines correspond (from top to bottom) to the evaluation of gi using the first one, two and 
three terms in the series eqn (20). Inset: the same quantity as a function of the scaled temperature r for 7 = 10. 



where 9 is the rapidity difference of the two particles. It is then easy to see that taking simultaneously 
the non-relativistic and weak-coupling limits of the sh-G model such that 
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0, c -> 00, g c = 4VX/H = fixed , (16) 



its iS-matrix (15) becomes identical to the ^-matrix (8) of the LL model [16, 28]. Notice that the coupling 
A does not need to be small, i.e. with this mapping one can study the LL model at arbitrarily large values 
of the dimcnsionlcss coupling 7. The mapping between the two models goes beyond the identity of their 
S- matrix: it extends both to their Lagrangians and TBA equations (details can be found in [28]). 

To perform the non-relativistic limit of the sh-G model, one has to express the real scalar field in the 
form 

Hx,t) = ^^(^(x,t)e-^ t + ^(x,t)e+^ t ) , (17) 

and, when the limit c — > 00 of the Lagrangian (or other expressions of (f>) is taken, of omitting all the 
oscillating terms [16, 28]. At equilibrium the expectation value of an operator O = 0(x) at temperature 
T and at finite density is given by 

[ ' Tr ( e -(ff-MA0/(fc B T)) ' 
In a rclativistic intcgrablc model the above quantity can be expressed as [29] 



<°> -e.s 7 (n £/<*>) 

n=0 _"L V=l / 



:V|o(o)i^) conn , (19) 



where / = 1/(1 + e e ^ e ^) and 9 = 61, . . . , 9 n (9 = 9„, . . . , #1) denote the asymptotic states entering 
the traces in (18). This formula contains both the pseudo-energy e{9) satisfying the TBA equations 
of the sh-G model and the connected diagonal form factor of the operator O. The latter is defined as 
(9 \0\ 9 )conn = J r (lim, )i ^o(0|0| 9 , 9 —iir + iV)) where rf = rj n , . . . , r)i and T in front of the expression 
means taking its finite part, i.e. omitting all the terms of the form r]i/r]j and l/^f where p is a positive 
integer. 

Since the form factors in the sh-G model are exactly known [30], one can compute local correlators in 
the LL model [16, 28]: 

<*' ***> = (?) "(£) " t a / (n £/<»>) ■■•■*>■ w 

Here f(p) = 1/(1 + e 5 ^-*) where e(p) is the solution of the non-relativistic TBA equations and 

<on»(fe}) - c ^ (^)"<o„n({^ = £}) (2D 
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FIG. 3: Two-body correlator gi as a function of 7 at T = for the STG gas. The solid line corresponds to the 
result obtained using the Hellmann-Feynman theorem (22), while the dashed lines (a), (b) and (c) are the results 
taking the first one, two and three terms in the series (20). The dotted lines (from the right) correspond to the 
asymptotic expansions (26) and (28). 



are the double limit (16) of the connected form factors. Note that this is a completely non-relativistic 
formula where only the form factors (21) have their origin in the sh-G model. However, following the 
arguments of [31] about the connection between the non-relativistic and relativistic form factors, the 
quantities (21) should also be derived solely from the LL model. 

Eqn (20) allows for the computation of the local correlators at equilibrium for the repulsive LL model, 
once the TBA equations with positive 7 are solved for the pseudo-energy e and the solution is inserted 
in (20) together with the form factor defined by (21). The way the method needs to be modified for 
negative 7 depends on the state for which we want to compute expectation values. If one is interested 
in computing the local correlators in the ground-state of the attractive LL model, one should use as a 
starting relativistic field theory the sine-Gordon model [32, 33] instead of the sh-G model. In the cold 
atom setup [9] this would correspond to a very slow, adiabatic switch from the positive to the negative 
side of the confinement-induced resonance such that the ground-state is reached. For a sudden switch 
when the meta-stable STG state is obtained, the Bethe roots remain real and their positions can be 
obtained from the Bethe ansatz equations of the LL model with negative 7. Since all the expectation 
values finally depend only on the Bethe roots, they can consistently be obtained by changing the sign of 
the coupling constant. Correspondingly our formula eqn (20) has to be used with negative 7. Of course, 
to be consistent one has to use the solutions of the TBA equations and the form factors with negative 7. 

Eqn (20) provides a series expansion in terms of the form factors. For small values of I7I we would have 
to sum many terms of the series. However, in the case of the STG gas we are interested in large values 
of I7I and as we discuss it in the next section, the series is rapidly converging and already the first few 
terms give very accurate results. 

In the next two sections we present the results for the STG local correlators and we compare them 
with available results from the literature. The two-body correlator (72 can be exactly determined [25, 26] 
via the Hellmann-Feynman theorem: 



For T > the theorem gives 



where the free energy F can be calculated from the TBA approach [see eqn (A4b)]. In dimensionless 
variables 

52(7, r) = ri ( a j g log(l + e~^) j , (24) 



where a = [ijk^,T with fi being the chemical potential. In [28] a compact form without derivatives was 
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FIG. 4: Three-body correlator g$ as a function of 7 at T = for the STG gas (we also plot for comparison the 
corresponding LL result). The STG solid line corresponds to the extension of [17] to negative values of 7, while 
the LL solid line is the plot of [17] for 7 > 0. The two dashed lines are the results taking the first one and two 
terms in the series (20). The dotted line corresponds to the large 7 asymptotic result (27). 



2 f j. 9(<l) „ 2 t [ dq ^ , 



derived 

^ = 7-y d «lf^-7 J iMl + e— )• (25) 

—00 —00 

As previously discussed, the above equations are also valid for the STG gas provided that 7 (and A and 
P in Appendix A) are chosen to be negative. 

For large 7 there are asymptotic expansions in the literature [16, 25]: e.g., for 172(7) at T = we have 
from (12) and the Hellmann-Feynman theorem 

9j ( 7 )^(l + ^\ W 
while for (73(7) we will use the leading order result of [25] which is invariant under the sign change 7 — > —7: 

16 7T 6 

157 



53(7) « -Ts- ( 27 ) 



Let us recall again our convention described at the end of Section II that for the STG results 7 denotes 
the absolute value of the effective coupling constant. 



IV. RESULTS FOR T = 



In this section we present our results for the two- and three-body correlators at zero temperature 
obtained by the method discussed in the previous section. To test the convergence of our series in Fig. 2 
we show the difference between our evaluation of 171(7) and the exact value 171 = (ip^ip) /n = 1. Although 
the convergence is slower than for the repulsive LL gas, for 7 = 10 the error using the first three terms 
in the series expansion is around 1%. 

In Fig. 3 we plot our results for 32(7) at T = 0. Our findings are compared with the exact result 
obtained using the Hellmann-Feynman theorem (22) and with the asymptotic expansions (26). We also 
plot the more accurate third order strong-coupling expansion obtained in [16]: 

»-£( 1+ 5 + ( M -§^). 

For intermediate values of 7 (in the region 5 < 7 < 25) our results are much better than the asymptotic 
results. A comparison with the repulsive LL results shows also that (72 for the STG gas in the intermediate 
region is significantly larger than the corresponding quantity for the LL gas. 

We show the results for 173(7) at T — in Fig. 4. Our results are plotted together with the result of [17] 
for the repulsive case. We also show the results for 5(3 obtained by extending the formula for 33 of [17] 
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FIG. 5: Two-body correlator for the STG gas as a function of the scaled temperature r for three different values 
of 7 = 10, 30, 100. The solid lines corresponds to the results obtained using the Hellmann-Feynman theorem (23), 
while the dashed lines are the results taking the first one, two and three terms of (20). 

to negative values of 7, corresponding to the STG gas. One can see that the first two terms of the series 
20 are in excellent agreement with the findings obtained from [17] already for 7 > 10. The comparison 
between the STG and the LL gas shows also that the former is much more subjected to three-body 
recombination: at 7 ~ 30, for example, g^ TG ^ /g^ L ^ ~ 3. 

V. RESULTS FOR T > 

In this section wc show our finite temperature results for gi and gz- To test the reliability of the series 
expansion (20) at finite temperature, we plot in the inset of Fig. 2 the deviations 1 — (71(7) from the exact 
result gi = 1 for a fixed value of 7. It is interesting to observe that for large coupling the convergence is 
faster than the LL result (not shown there) . From the figure we can also see that the convergence slightly 
improves with increasing temperature. 

In Fig. 5 #2 is plotted as a function of the scaled temperature for three different values of 7, together 
with the results obtained from the Hellmann-Feynman theorem (23) which agree very well with our 
expansion. It can be seen in Fig. 5 that g2 is approximately constant for temperatures T < Td: e.g. for 
7 = 100 wc find gi(r = V)j gi(r = 0) « 1.1. For higher values of 7 the temperature effects get relevant at 
slightly lower temperatures. Comparing g^ (r) for the STG and LL gases it is also possible to see that for 
the STG gas temperature effects become important at slightly higher temperatures than for the LL gas. 

In Figs. 6-7 we finally plot our results for g 3 for the STG gas at finite temperature. In Fig. 6 we show 
(73 for the STG gas as a function of 7 for two different scaled temperatures r = 1 and r = 10. In the 
inset we plot 33 for the STG and LL gases as a function of r for 7 = 10 and 7 = 30: it is clear that for 
the STG gas 33 is larger than the corresponding LL value and that in both cases temperature effects arc 
present at T > Td (for the same 7 significant deviations from the T = results start at slightly higher 
temperatures for the STG gas). 

In Fig. 7 we plot the ratio g^ / g^ between the three-body correlators of the STG and LL gases 

as functions of 7 for three different scaled temperatures. In the inset of the same figure g^ TG ^ / g^ L ^ is 
plotted for two values of 7 as a function of r. From Fig. 7 one can see that for T = Td (t = 1) this ratio 
is very large for intermediate values of the coupling and it is not very different from the T = result. For 
large values of 7 the ratio becomes smaller: at 7 = 30 the ratio is ~ 3 for t = 1 and ~ 1.5 for t = 10. 

VI. CONCLUSIONS 

We studied the local correlations in the super Tonks-Girardcau gas, a highly excited, strongly correlated 
state obtained in quasi onc-dimcnsional Bosc gases when the scattering length is tuned to large negative 
values using a confinement-induced resonance. After introducing the Lieb-Liniger model we discussed 
the main properties of the super Tonks-Girardeau gas which was recently realized in the experiment 
reported in [9]. We focused on the computation of the local correlators: using a relation with a relativistic 
field theory, we obtained results for the two-body and three-body local correlators at zero and finite 
temperature. At T = we showed that the three-body correlator agrees with the extension of the results 
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FIG. 6: Three-body correlator go, for the STG gas as a function of 7 for two different scaled temperatures r = 1, 10. 
The two dashed lines for each value of r are the results taking the first one (top) and two (bottom) terms in the 
series (20). Inset: gs for the STG and LL gases as functions of the scaled temperature r for two different values 
of 7 = 10,30. 
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FIG. 7: Ratio between the three-body correlator gs, of the STG gas and the corresponding quantity for the LL gas 
at three different scaled temperatures r = 0, 1, 10 (from top to bottom). The r = solid line is obtained using 
the results of [17]. Inset: the same quantity as a function of r for two different couplings 7 = 10, 30. In both cases 
the dashed lines are obtained using two terms of the series (20). 




of Cheianov et al. [17] obtained for the ground-state of the repulsive Lieb-Liniger gas, to the super 
Tonks-Girardeau state. At finite temperature, the three-body correlator for intermediate values of the 
coupling constant has a very weak dependence on the temperature up to the degeneracy temperature Td- 
We also showed that the value of g% for larger temperatures at even relatively small coupling constants 
7 is rather similar to the corresponding value of the repulsive Lieb-Liniger gas. 
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Appendix A: Thermodynamical Bethe ansatz equations for the Lieb-Liniger model 

In this Appendix we summarize the TBA equations describing the equilibrium properties of a repulsive 
LL gas. In the limit N — > 00, L — > 00 with the density n fixed, the discrete energy levels of the system 
get encoded in an energy level density function p(p) and in the density p^(p) of the occupied levels. The 
ratio between the two densities p and p^ 1 ' defines the pseudo-energy e(p) through the relation 

p{p) = 1 + e^ , (Al) 



and this quantity, together with the densities, satisfies a coupled set of integral equations. Using the 
rescaled quantities 

g= — , a=— — , g{q) = hp(nhq) , (A2) 
nn fce ± 
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where [i is the chemical potential, T is the temperature and fee is the Boltzmann constant, the equations 
read 

oo 

q 2 f dq' 27 



and 



— OO 

CXD 

9 ( g ) = _L+ / d q'- ^ ^W) (A3b) 

yW 2tt 7 1 {q-q'f+l 2 l + e £ («') ' V ' 



-dg. (A3c) 



1 + e £ w 



The physical parameters of the problem are A, T and n, but only the dimensionless combinations 7 and 
t enter the results. The chemical potential (or the dimensionless fugacity-like parameter a) gets fixed by 
the constraint given by the last equation. 

Once the TBA integral equations (A3) are solved, the ground-state energy E and the free energy F of 
the system are expressed as 

00 00 
E ~ fdpf- pW (p) = ^n* [ dq -^q 2 , (A4a) 



L / 2m v ' 2m / 1 + e £ (?) 



= nk B T la--^ J dq log (l + e -eC «>) J . (A4b) 



At zero temperature the energy level density gets a compact support, i.e. it is different from zero only 
on an interval (which we denote by [-B, B]) and, correspondingly, the TBA equations simplify. Applying 
a different rescaling, 

* = = ( A5) 

we arrive at the LL integral equations 

1 

1 = 1 J v {k) dk , (A6a) 



-1 



1 



-1 

while the ground-state energy is given by 

B 

„2 fc2 



1 f dk' 2/3 , lM 

"W = ^ + / 27 (t -tr + ^ ^»' (A6b » 



f ^ / " £" 3 (?) / dMk)k * s ss" 3eW ■ (A7) 

—23 -1 
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